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Abstract 

We present new results of a deep optical imaging survey using a narrow band 
filter (A^_B921) centered at A = 9196 A together with B, V, R, i', and z' broadband 
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filters in the sky area of tlie Subaru Deep Field, which has been promoted as one of 
legacy programs of the 8.2 m Subaru Telescope. We obtained a photometric sample 
of 58 Lya emitter candidates at z ^ 6.5 - 6.6 among ~ 180 strong iVi?921-excess 
{z' — NB921 > 1.0) objects together with a color criterion of i' — z' > 1.3. We then 
obtained optical spectra of 20 objects in our A^i?921-excess sample, and identified 
at least nine Lya emitters at 2 ~ 6.5 - 6.6, including the two emitters reported by 
Kodaira et al. (2003, PASJ, 55, L17). Since our Lya-emitter candidates are free from 
strong amplification of gravitational lensing, we are able to discuss their observational 
properties from a statistical point of view. Based on these new results, we obtained 
a lower limit of the star- formation rate density of Psfr — 5.7 x 10~^ /io.7 Mq yr~^ 
Mpc~^ a.t z^ 6.6, being consistent with our previous estimate. We discuss the nature 
of star-formation activity in galaxies beyond z = 6. 

Key words: cosmology: observations — early universe — galaxies: formation 
— galaxies: evolution — galaxies: starburst 

1. INTRODUCTION 

A simple method to understand the formation process of galaxies is to find a sample 
of very young galaxies at high redshift, and then to investigate their observational properties 
in detail. Theoretical models based on the hierarchical clustering scenario suggest that first- 
generation (i.e.. Population III) objects could be born around 2; ~ 30 (~ 0.5 Gyr after the big 
bang), and then galactic systems with masses higher than IO^^Mq could be assembled after 2; ~ 5 
- 10 (~ 1 Gyr after the big bang) although their co-moving number density could be significantly 
smaller than that in the present day (e.g., Ostriker, Gnedin 1996). It is, thus, important to 
search for star-forming galaxies beyond redshift 5 to probe the cosmic star-formation history 
from very high redshift to the present day. 

Surveys for such high- 2; galaxies have been made mostly by the optical color-selection 
technique [e.g., Steidel et al. (1996), Madau et al. (1996) see for the recent discovery of a 
galaxy at z = 10, Pello et al. (2004a); note that the detection of Lya emission line has been 
still in debate; see Weatherley et al. (2004); Pello et al. (2004b)]. Indeed, recent wide- field, 
deep imaging surveys using the Advanced Camera for Surveys on board the Hubble Space 
Telescope surely found a sample of galaxies at 2 ~ 6 (Giavalisco et al. 2004b; Dickinson et 
al. 2004; Bouwens et al. 2004; Stanway et al. 2004a, 2004b). In addition, recent advances 
in deep optical imaging capability with 8-10 m class telescopes have enabled new searches for 

^ Based on data collected at the Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 

^ This work has been done with a collaboration between the Subaru Deep Field Project led by the association 
of builders of the Subaru telescope and the common-use, Intensive Program (S02A-IP-2) led by Y. Taniguchi. 
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star- forming galaxies beyond redshift 5 (see for review, Taniguchi et al. 2003b; Spinrad 2003). 
In particular, imaging surveys using narrow-passband filters have proved to be an efficient way 
to find such galaxies (e.g., Hu et al. 2002, 2004; Kodaira et al. 2003; Maier et al. 2003; Rhoads 
et al. 2003, 2004). These surveys can probe objects with a high star- formation rate, even if 
they are too faint to be detected in color-selection procedures. Therefore, these two methods 
cooperate with each other to investigate young, star- forming galaxies at high redshift. It is also 
worthwhile noting that recent imaging-spectroscopic surveys with a grism are also capable of 
finding high-2; LAEs. A recent such trial succeeded to find a LAE aX z = 6.518 (Kurk et al. 
2004). A serendipitous discovery of a LAE at z = 6.545 was also reported (Stern et al. 2004). 
For a review of discovery methods, see Taniguchi et al. (2003b). 

The Subaru Telescope team officially started a large-scale deep-survey program in 2002 
April; the Subaru Deep Field (SDF) project. Several pilot papers related to the SDF project 
have already been published (Maihara et al. 2001; Totani et al. 2001a, 2001b, 2001c; Ouchi et 
al. 2001, 2003, 2004a, 2004b; Shimasaku et al. 2003; Kashikawa et al. 2003). The main aim 
of this project is to investigate the formation and evolution of galaxies from very high redshift 
[z ~ 7) to the present day {z ~ 0) based on a set of very deep optical and near infrared (NIR) 
imaging data and follow-up spectroscopy in both the optical and NIR. In addition, other multi- 
wavelength observations will be added in the future. The overall design of the SDF project will 
be given elsewhere (Kashikawa et al. 2004a). 

Since the SDF data set is very deep in both the broad and narrow filter bands (see 
table 1), we will be able to contribute to the progress in searches for such high-z galaxies based 
on the above two methods. Using our SDF data taken in 2002, we have already found two 
Lya emitters (LAEs) at 2 = 6.578 and z = 6.541 (Kodaira et al. 2003). The former object is 
the most distant Lya emitter found in optical surveys, SDF J132418. 3+271455; note that the 
second most distant known is HCM 6A at z = 6.56 (Hu et al. 2002). In this paper, we present 
a summary of our search for LAEs at z ^ 6.6 made during these two years. We adopt a flat 
universe with ^ma,ttev = 0.3, = 0.7, and h^ j = Hq/ (70 km s~^ Mpc~^) throughout this paper. 
We use the AB system for optical magnitudes. 

2. OBSERVATIONS 

2.1. Optical Deep Imaging 

We have carried out a very deep optical imaging survey in the SDF centered at a(J2000) 
= 13^ 24'^ 38^9 and 5(J2000) = +27° 29' 25".9 (Kashikawa et al. 2004a). Optical imaging 
was made in the B, V, R, i', z', NB816, and A^i?921 bands on a central 34' x 27' area of the 
SDF with Suprime-Cam, which consisted of 5 x 2 CCDs of 2k x 4k pixels, with a pixel scale 
of 0."202 pixel"^ (Miyazaki et al. 2002) on the 8.2 m Subaru Telescope (Kaifu et al. 2000; lye 
et al. 2004). The data were collected in several observing runs during a period between 2001 
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and 2003. A summary of the imaging observations is given in table 1. 

Two narrow band filters were used to find LAEs at z ^ 5.7, and z ~ 6.6. Other narrow 
band filters, NB70A and NB711, were also used in the SDF to map the SDF. Early results 
with the use of NB711 were reported by Ouchi et al. (2003) and Shimasaku et al. (2003), 
and those with the use of NB7QA were by Shimasaku et al. (2004). In this paper, we present 
results obtained with the use of NB921 filter centered on Ac = 9196 A with a passband of 
AA(FWHM) = 132 A, together with broad-band data. The central wavelength corresponds to 
a redshift of 6.56 for Lya emission. Details of the filter system used in the SDF project are 
given in Kashikawa et al. (2004a). The central wavelength corresponds to a redshift of 6.56 for 
Lya emission. Results obtained by using our NB816 data will be reported elsewhere. 

The individual CCD data were reduced and combined using both our own data-reduction 
software (Yagi et al. 2002) and IRAF. The combined images for individual bands were aligned 
and smoothed with Gaussian kernels to match their seeing sizes. The PSF FWHM of the final 
images is 0."98. The exposure times and limiting magnitudes are listed in table 1. Photometric 
calibrations were made using the usual spectrophotometric standard stars. Source detection 
and photometry were performed using SExtractor (Bertin, Arnouts 1996) version 2.1.6. The 
iVi?921-band image was chosen to detect candidate objects, and we limit the object catalog to 
NB921 < 26.54, the 3(T limiting magnitude. For each object detected in the NB921 image, the 
z', and NB921 magnitudes were measured on a common aperture of 2."0 diameter. In total, 
~ 82000 objects were detected down to iV5921 = 26.54. 

The bandpass of the NB921 filter is completely included in that of the z' filter. Therefore, 
when the Lyman break (Arest = 912 A), or nearly zero continuum flux at Arest(Lya) = 1216 A is 
redshifted to the z' window, Lyman-break galaxies could show some excess in the z' — NB921 
color, even if they have little Lya emission. In order to examine this possibility, in figure 1, we 
show the z' — NB921 color as a function of the redshift for a galaxy without Lya emission; the 
model galaxy is generated by using GALAXEV (Bruzual, Chariot 2003) with parameters of 
Z = 0.02, r = 1 Gyr, and the age = 1 Gyr. Even for such a Lyman break galaxy without Lya 
emission, the z' — NB921 color becomes as high as 0.9 mag at 2; ~ 6.4. Therefore, we required 
that the z' — NB921 color exceeds 1.0 mag for strong Lya emitters at z ~ 6.6; i.e., 

z' -NB921>1.0. (1) 

Together with this criterion, we also used the following two more selection criteria to select 
A^i?921-excess objects: 

z' - NB921 > 3a (2) 

and 

NB921<26.0{5a). (3) 
Although our 3a detection limit for NB921 was 26.54, we adopted the above 5a limit to secure 
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the selection of A^i?921-excess objects. Using these criteria, we obtain a sample of 185 NB921- 
excess objects. In figure 2, we show a diagram between z' — NB921 and A^i?921, where these 
criteria are also shown. 

In order to reduce contamination from foreground objects that are free from absorption 
by the intergalactic medium, we adopted another color criterion, 

i'-z >1.3, (4) 

together with that LAE candidates are not detected in S, V, and R band images (less than 
3cr in each band). In order to show this color criterion, we plot a diagram between z' — NB921 
and i' — z' in figure 3. We also show the color evolution of a model galaxy as a function of the 
redshift; colors of a starburst galaxy were calculated by using the population synthesis model 
GALAXEV (Bruzual, Chariot 2003) where we adopted the r = 1 Gyr model with an age of 
t = 1 Gyr and adding some emission lines, e.g., Lya, [O ll], and [O ill]. The emission-line 
luminosities were calculated in the following way: (i) the H/5 luminosity was calculated from 
the ionizing photon production rate, A^Lyc photons s~^, using the relation 

L(H/?)(erg s'^) = 4.76 x lO'^^iVLycls"'), (5) 

and (ii) the other emission-line luminosities are calculated using the relative luminosity to H/3 
luminosity tabulated in PEGASE (Fioc, Rocca-Volmerange 1997). 

In order not to miss possible faint LAEs, we included all objects with z' — NB921 > 1.0 
even fainter than i' = 27.87 (~ 2a limiting mag) in our LAE sample. Part of them nominally 
had i' — z' > 1.3, although we did not use the i' — z' color criterion for objects fainter than 
z' ~ 26.57 (= 27.87 — 1.3). By using all of these criteria, we obtained a photometric sample of 
58 LAE candidates. Their basic properties are given in table 2. In this table, we also give a pure 
z'-band magnitude, 2:9500; for each object. This magnitude was corrected for the contribution 
of the Lya flux, estimated by using the NB921 magnitude. The corresponding pure ^'-band 
flux was estimated as 

/9500 = i^K'fiz') - 0.7AAjvB92i/(iV5921)]/AA^/, (6) 

where f{z') and f{NB921) are the observed z' and A^i?921 fluxes, respectively, AA2' is the 
bandpass of the z' filter in units of A (= 960 A), AXj\fB92i is that of the NB921 filter (= 132 
A), AA^? is the effective bandpass of the z' filter for Aobs > 9260 A (=270 A), and the numerical 
factor of 0.7 is the relative transmittance of the NB921 filter with respect to the z' filter. 

The Lya flux corrected for the contribution of UV continuum emission at wavelengths 
longer than 1216 A for each object was estimated by using 

/Lya,image = /(A^5921) - /a(2950o) X AA^B92l/2. (7) 

Here, we note that this sample was slightly different from that discussed in Kodaira et 
al. (2003) and Taniguchi (2003), because the photometric data used in their papers were based 
on data obtained before 2002 May and their LAE selection criteria were also different from 
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those used in this paper. The previous LAE selection criteria adopted in Kodaira et ah (2003) 
were: (1) z'-NB921 > 1.0 and2'-^'> 1.3 for objects with i' > 28.0, and (2) z'-NB921 > 1.0 
for objects with i' < 28.0. Note that Taniguchi (2003) adopted z' - NB921 > 0.9 because 
the photometric sample at that occasion was tentatively selected for our follow-up optical 
spectroscopy. Note that the photometric catalog used in this paper is the final one described 
in Kashikawa et al. (2004a) and thus the selection procedure of A^i?921-excess objects is also 
our final one. 

The observed equivalent width of Lya emission for each galaxy in our photometric 
sample is given in table 2. Our NB921-excess criterion given in (1) nominally corresponds 
to the cutoff for the Lya-emission equivalent width, EWohsO^Jct) > 200 A. However, since the 
actual equivalent width of Lya emission should be obtained from EWohs(J-'JC() = /Lya,imago/ /gsoo, 
the NB921-excess criterion given in (1) gives the equivalent width cutoff of EWohsO^Jc^) > 50 
A, corresponding to the rest-frame equivalent width, EWo(Lya) > 7 A. However, most LAE 
candidates in our sample tend to be faint in the z' band. Further, in order to secure the 
selection of LAE candidates, we also adopted the criterion (2): z' — NB921 > 3a. Therefore, 
LAEs with much stronger Lya emission could be detected in our analysis; our nine LAEs have 
EWobsC^-ja) > 130 A (see table 2). 

The effective area used to search for A^i?921-excess objects was 875.4 arcmin^. The 
FWHM half-power points of the filter corresponded to a co-moving depth along the line of 
sight of 40.9 Hq I Mpc {z-^i^ ^ 6.508 and Zmax ~ 6.617; note that the transmission curve of our 
NB921 filter has a Gaussian-like shape). Therefore, a total volume of 217200 h^^ Mpc^ was 
probed in our A^5921 image. 

2.2. Optical Spectroscopy 

In order to investigate the nature of LAE candidates found in our optical-imaging survey, 
we obtained optical spectra of 20 objects in our photometric sample of LAEs; note that the nine 
galaxies for which optical spectroscopy was made were randomly selected from a photometric 
sample. We used the Subaru Faint Object Camera And Spectrograph (FOCAS; Kashikawa et 
al. 2002) on 2002 June and 2003 May and June. A journal of our spectroscopy is summarized 
in table 3. Details of this spectroscopy will be given elsewhere. 

Our spectroscopy was made with the multi-object slit (MOS) mode using the following 
grating sets: (i) a 300 lines mm"^ grating with an 058 order cut filter, and (ii) an Echelle with 
a z' filter. An 0."8-wide slit was used in both settings The wavelength coverage was ~ 6000 
A to 10000 A for case (i), while it was ~ 8000 A to 10000 A for case (ii). The spectroscopic 
resolution was 9.0 A at 9200 A, i? ~ 1020, for case (i) while 6.3 A at 9000 A, i? ~ 1430, for case 
(ii). The spatial sampling rate was 0."3 pixel"^ after 3-pixel, on-chip binning in both cases. 
The flux calibration was made with spectra of spectroscopic standard stars, Hz 44 and Feige 
34. 
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In figure 4, we show the distributions of the NB921 magnitudes for the photomet- 
ric sample (58 objects) and the spectroscopic sample (20 objects); note that those of the 
spectroscopically-confirmed LAE sample (9 objects), given in subsection 3.1, are also shown. 
In figure 5, we show the spatial distributions of the above three samples. Details on this issue 
will be discussed by Kashikawa et al. (2004b, in preparation). 

3. RESULTS 

3.1. Spectroscopic Identification 

Our photometric sample of LAE candidates consisted of 58 objects. We obtained optical 
spectra of 20 objects among them. We detected single emission lines in 14 objects, but did not 
detected any emission-line feature in the remaining six objects, although one shows a marginal 
continuum feature. 

The fourteen emission-line objects with a single emission line at A ~ 9200 A may be 
either an LAE a.t z ^ 6.6 or an [O ll] emitter at z ~ 1.47 (e.g.. Stern et al. 2000). LAEs at 
high redshift show a sharp cutoff at wavelengths shortward of the line peak because of H I 
absorption by gas clouds in the system and intergalactic H I gas (e.g., Hu et al. 2002; Ajiki et 
al. 2002). Another important spectral feature of LAEs is that there is little continuum emission 
shortward of the Lya line (Hu et al. 2002; Kodaira et al. 2003). Such a continuum break is 
actually seen in our two objects, SDF J132415.7+273058 and SDF J132418.3+271455, which 
are already identified as LAEs at 2; ~ 6.6 (Kodaira et al. 2003). However, the other objects show 
no continuum feature at wavelengths longer than the line peak, and thus we cannot identify 
LAEs solely by using the continuum-break feature. 

Among the remaining emission-line objects, the following two have a narrow line width: 
6.5 A for SDF J132518.8+273043, and 5.5 A for SDF J132522.3+273520. If the emission line 
were [O ll] emission, the redshift would be 2; ~ 1.47. Since the [O ll] feature is a doublet line of 
[O ll]A3726.0 and [O ll]A3728.8, the line separation would be wider than 6.9 A at 2 = 1.47 and 
the two components could be resolved in our Echelle spectroscopy. Further, if the line were H/3, 
[O lll]A4959, [O lll]A5007, or Ha line, we would detect some other emission lines in our spectra. 
Therefore, these two narrow-line objects must be a LAE (e.g., Taniguchi et al. 2003a). 

As for the other ten emission-line objects, we found that the following five objects show 
an asymmetric profile with a sharp cutoff at wavelengths shortward of the line peak: SDF 
J132352. 7+271622, SDF J132353. 1+271631, SDF J132408.3+271543, SDF J132418.4+273345, 
and SDF J132432. 5+271647. Therefore, the nine objects discussed above appear to be nice 
candidates for LAEs at z ^ 6.5 - 6.6. The remaining five emission-line objects show either a 
symmetric emission-line profile or poor-S/N spectra. 

In order to secure our spectral classification, we estimated the flux ratio between /red 
and /blue, where /red is the flux at wavelengths longer than the emission- line peak, while fuue 
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is that at wavelengths longer than the emission- line peak. In table 4, we give our results. In 
figure 6, we show the frequency distribution of log /rcd//biuc for the sample of 14 emission-line 
objects. In this figure, we find bi-modal distributions in log /rcd//biuc; i-e., the nine objects 
thought to be nice LAE candidates have /rod//biuc > 1, while the remaining five objects have 
/rcd//biuc < 1- Some other observational tests have been proposed to identify LAEs at high 
redshift (e.g., Rhoads et al. 2003; Kurk et al. 2004); the S/N in our optical spectra prevents 
us from applying such tests. 

In total, we identified nine LAEs whose emission-line shapes show a sharp cutoff at 
wavelengths shortward of the line peak. They are listed up in the upper part of table 2. Their 
optical spectra are shown in figure 7. 

Since the remaining five emission-line objects show either a symmetric emission-line 
profile, or poor-S/N spectra, we did not include them as our final LAE sample, although we 
cannot rule out the possibility that they are LAEs at 2; ~ 6.6. They are listed in the middle 
part of table 2, with a Lya redshift, if they are a LAE. Their spectra are shown in figure 8. 

We found no emission-line feature in the six candidates in our spectroscopic sample (Nos. 
10, 11, 13, 14, 15, and 20 in table 2). The estimated noises in our FOCAS spectra in these six 
objects range from 2 x 10"^^ erg s~^ cm^ to 3 x 10~^^ erg s~^ cm^ A~^. Taking account of 
these noises, we could detect their Lya fiuxes, expected from our NB921 imaging, if they had 
an unresolved Lya emission line in our FOCAS Echelle spectra (6.4 A resolution), because 
their peak fiuxes would range from 3 x 10^^^ erg s~^ cm^ A~^ to 6 x 10~^^ erg s~^ cm^ A~^. 
If the line had some broad feature, or if the line was recorded in some noisy spectral parts, it 
seems that it would have been quite difficult to detect unambiguous emission-line feature for 
these six objects in our FOCAS spectroscopy. Since there was no emission-line feature in these 
six objects, we cannot judge whether they are LAEs or low- 2; emission-line sources. Another 
possibility is that they might be LBGs with weak Lya emission. In any case, we need future 
very deep optical spectroscopy for these objects. 

Thumbnail images of the nine LAEs are shown in figure 9; their photometric properties 
are given in table 2. The redshift given in the last column of table 2 was estimated from the 
peak of the Lya emission line. The measurement error of redshift was estimated to be ± 0.002. 
The observed fiux of Lya emission of each object is given in table 5. The Lya line width is also 
given in this table. The line width ranges from 180 km s~^ to 460 km s~^, being comparable to 
those of previously identified LAEs beyond z = 5. Our spectra show no evidence of N v A 1240 
emission for all the nine LAEs. These properties suggest that they are star-forming galaxies, 
rather than active galactic nuclei. LAEs with FWHM > 300 km s~^ (four among the nine) may 
experience a superwind in some cases (e.g., Dawson et al. 2001; Ajiki et al. 2002). 
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3.2. A Fraction of LAEs 

Our optical spectroscopy has confirmed that 9 among 20 candidates are real LAEs at 
z ~ 6.6. Five among the remaining 11 objects are single-line emitters with a symmetric emission- 
line profile. They may be either a LAE at 2; ~ 6.6 or an [O ll] A 3727 emitter at 2; ~ 1.47. 

Since it is safe to keep the five single-line objects as unclassified objects, we obtain from 
our spectroscopy that the fraction of reliable LAEs in our photometric sample is /(LAE) = 
9/20 = 0.45. This can be regarded as a lower limit in our study. If we include nominally 5 
single-line emitters into the LAE sample, we obtain /(LAE) = 14/20 = 0.70. The real fraction 
of LAEs in our study seems to be in the range of /(LAE) ~ 0.45 - 0.70. 

3.3. Morphological Structures of LAEs 

Since our final PSF size is 0."98, it seems difficult to investigate the structural properties 
of the nine LAEs. However, some NB921 CCD frames were obtained under 0."5 - 0."7 conditions 
during our observing runs. Using only these data, we made a high-resolution NB921 image with 
a PSF of 0."71. The integration time for this image was 37215.7 seconds, ~ 70% of the total 
integration time (see table 1). 

In figure 10, we show newly made NB921 images of the nine LAEs. Their azimuthal 
radial profiles are also shown in figure 11. Their sizes (FWHM) range from 0.81 arcsec to 1.02 
arcsec, as given in table 6. Although we do not conclude that all LAEs are spatially resolved 
in our NB921 images, some of them (e.g., Nos. 4 and 8) appear to be spatially extended. 
After deconvolution with the PSF size, we obtain their sizes between 0.39 and 0.73 arcsec, 
corresponding to ~ (2 - 4) h^^ kpc. It is thus suggested that gaseous matter around LAE 
host galaxies are spatially extended up to several kpc at 2; ~ 6.6 in some cases. Imaging with 
Advanced Camera for Surveys on board the Hubble Space Telescope with the 2;(F850LP) filter 
will be important to investigate their detailed Lya morphologies. 

3.4. Comments on Possible Amplification by Gravitational Lensing 

The SDF is a so-called blank field and thus there is no apparent cluster of galaxies 
known to date at low and intermediate redshifts in our field. In this respect, our survey may 
not suffer from strong amplification by gravitational lensing, unlike several surveys made with 
the help of gravitational lensing (Ellis et al. 2001; Hu et al. 2002; Santos et al. 2004; Kneib 
et al. 2004; Pello et al. 2004a). However, any high- 2; objects could suffer from gravitational 
lensing because of a larger optical depth for lensing (Wyithe, Loeb 2002; Shioya et al. 2002). 
Therefore, it is possible that some LAEs found in our survey could be gravitationally amplified 
by a foreground galaxy lying close to their lines of sight. 

In order to check such possibilities, we examined carefully our i' and z' images for all 
the nine objects. We also examined our B, V, and R images to check whether or not there 
are any foreground galaxies. As shown in figures 9 and 10, there are some galaxies that could 
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be foreground galaxies around each LAE. For example, a probable foreground object is found 
to be located at 0."8 NW of SDF J132418.3+271455. We listed all such foreground galaxies 
within a radius of 4". We then evaluated the magnification factor for each LAE in the following 
way (Shioya et al. 2002; see also Yamada et al. 2003): (1) We estimate its photometric redshift 
and rest-frame -B-band luminosity based on our optical broad-band photometric data for each 
foreground galaxy, if any. (2) We estimated its stellar velocity dispersion using the TuUy-Fisher 
relation and the -B-band luminosity estimated above. (3) We estimated the magnification factor 
based on the singular isothermal sphere model for gravitational lensing. We found that the 
amplification factor by gravitational lensing is smaller than 1.1 for all the cases. It is also 
noted that there is no counter image for each LAE within the limiting magnitude, and thus 
the gravitational amplification factor should be less than a factor of ~ 2, being consistent with 
the above factor. Therefore, we conclude that our LAE sample did not suffer from strong 
amplification by the gravitational lensing. This allows us to perform simple statistical analyses 
for our sample of LAEs, such as the number density and the star-formation rate density at 
z = 6.6. 

4. DISCUSSION 

4.I. The Number Density of LAEs at z ^ 6.6 

First, we estimated the number density of LAEs at 2; ~ 6.6 and compared it with those 
at other redshifts. Our sample of LAE candidates consists of 58 objects. Then, our optical 
spectroscopy found that nine objects are LAEs at 2 ^ 6.5 - 6.6 among the 20 objects in our 
spectroscopic sample. Therefore, from a statistical point of view, ^ 45% of our photometrically 
selected candidates are expected to be LAEs at 2 ~ 6.6; note again that the nine galaxies for 
which optical spectroscopy was made were randomly selected from the photometric sample. The 
limiting magnitude of our sample, NB921 = 26.54, corresponds roughly to /(Lya) ~ 4.1 x 10~^^ 
erg s~^ cm~^, or L(Lya;) ~ 2.0 x 10^^ h^j erg s~^ for our adopted cosmology. Since the survey 
volume of our survey is 217200 h^^ Mpc'^, we obtain the number density of LAEs brighter than 
this limiting luminosity to be n(LAE) ~ 58 x 0.45/217200 ~ 1.2 x 10"^ h^ ^ Mpc'^ 

It is interesting to compare our result with previous studies on LAEs at high redshift. 
Note, however, that the detection completeness is not taken into account in later discussion 
because it is difficult to estimate it for each survey. Ouchi et al. (2003) have made a wide-field 
survey of LAEs at 2; ^ 4.9. From their data, we estimate that the number density of LAEs 
at z ~ 4.9 brighter than L(Lya) = 2.4 x 10^^ /^-2 g-i jg ^ 2 x 10""^ 7 Mpc'^ (M. Ouchi, 
private communication), which is close to that for our z ~ 6.6 LAEs. Rhoads and Malhotra 
(2001) also carried out a similar survey for z = 5.7 LAEs; the limiting luminosity of their survey, 
!=s 5.3 X 10^^ erg s^^, is brighter than ours (see also Rhoads et al. 2003) . They found in 
their sample the number density of LAEs to be ~ 4 x 10"^ 7 Mpc~^. If we limit our z ~ 6.6 
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LAE sample to Rhoads et al.'s (2003) limiting luminosity, we obtain ~ 2 x 10~^ /ig 7 Mpc~^. 
Similarly, we obtain ~ 4 x 10^^ h^j Mpc~^ for Ouchi et al.'s (2003) sample. These results 
indicate that the number density of LAEs does not change significantly from z ~ 4.9 to 6.6. 
Combined with a result given in Ouchi et al. (2003) that the number density of 2; ~ 4.9 LAEs is 
not clearly different from that of 2 = 3.4 LAEs in Cowie and Hu's (1998) sample, we conclude 
that LAEs do not evolve in number from z ~ 3.4 to 6.6 within the uncertainties in the data. 

4-2. The Star- Formation Rate of LAEs at z ^ 6.6 

4-2.1. The Star-Formation Rate Estimated from Lya Luminosity 

Next, we estimate the star-formation rate of the LAEs at z ^ 6.6. We have two kinds 
of SFRs: one is based on our optical spectroscopy and the other is based on our continuum- 
subtracted NB921 data. 

First, we discuss the SFRs based on our spectroscopy. The observed Lya flux, Lya 
luminosity, and star- formation rate, SFR(Lya), of each LAE are summarized in table 5. Note 
that the SFR(Lya) is estimated by using the relation (Kennicutt 1998; Brocklehurst 1971) 

SFRiLya) = 9.1 x 10"^^L(Lya) Mq yr-\ (8) 

where the Salpeter initial mass function with (miowor, '"^upper) = (0.1 Mq, 100 Mq) is adopted. 
The SFRs obtained for the nine LAEs range from ^ 3 Hqj Mq yr~^ to ~ 9 h^^ Mq yr~^ with 
an average of 5.7 ± 2.3 h^^ Mq yr~^. These values are comparable to those of LAEs at 2; ~ 
5.1 - 5.8 (e.g., Taniguchi et al. 2003b and references therein). The total spectroscopic SFR for 
the nine LAEs amounts to 

^^^spec^Lya) ~ 49.4 Hq^j Mq yi-\ (9) 

It should be mentioned that the SFRs estimated above are lower limits because it is quite likely 
that a blue half or more of the Lya emission may be absorbed by dust grains in the galaxy, 
itself, and by the intergalactic H I gas (e.g., Hu et al. 2002; Haiman 2002). Therefore, it would 
be desirable to study the rest-frame UV-optical continuum (Hu et al. 2002; Kodaira et al. 
2003). 

Second, we discuss SFRs based on the continuum-subtracted NB921 data; 
SFR^^^^^{Lya). Our results are summarized in table 7; note that we give the photometric SFRs 
for all 58 LAE candidates. The frequency distributions of Lya luminosity for all 58 LAE can- 
didates are shown in figure 12. The Lya luminosity of the two sources, SDF J132353. 1+271631 
(No. 2) and SDF J132518. 8+273043, is given as an upper-limit. This is due to the fact that 
the bandpass-corrected z' flux (/9500; see table 2) is estimated to be a bit higher because of low 
signal-to-noise ratios. Although these two upper limit data are present, the SFRs of the nine 
LAEs confirmed by our spectroscopy range between 3 and 10 /i(^7 MQyr"^, being similar to 
that obtained in our spectroscopy. The total photometric SFR for the nine LAEs is estimated 
to be 
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^^^image^Lya) ~ 58.5 /la? Mq yr-\ (10) 

being nearly the same as that obtained by the spectroscopy. In figure 13, we compare the 
observed fluxes of Lya emission obtained from spectroscopy and imaging for the nine LAEs. 

4-2.2. The Star Formation Rate Estimated from UV Luminosity 

We now estimate another SFR derived from the UV continuum luminosity for our sample. 
The observed z' magnitude can be converted to a UV continuum luminosity at A = 1260 A. 
Using the relation (Kennicutt 1998; see also Madau et al. 1998), 

SFR(UV) = 1.4 X 10^2^L^(UV) Mq yr"^ (11) 

where L;^(UV) is the UV continuum luminosity in units of erg s~^ Hz^^, we estimate the SFR 
based on the rest-frame UV (A = 1260A) continuum luminosity for each object. The results are 
summarized in table 7. Note that the same initial mass function for SFR{Lja) was adopted 
in this estimate. The frequency distributions of UV continuum luminosity for all 58 LAE 
candidates are shown in figure 14. 

4.2.3. Comparisons between SFR(Lya) and SFR(UV) 

We then compared the SFRs derived from the Lya and the UV photometric data. The 
result is shown in figure 15; also see table 7. As expected from previous studies (e.g., Hu et al. 
2003; Kodaira et al. 2003; Ajiki et al. 2003), we find that SFR{\]Y) tends to be higher by a 
factor of 5, on average, than S'-Fi?'™'^s*^(Lya). The obtained factor seems to be higher than the 
previous results; e.g., by a factor of 2. 

4.3. The Star Formation Rate Density of LAEs at z 6.6 

We now discuss the star-formation history beyond redshift 6 based on our data. We can 
estimate the total star-formation rate of 58 LAEs in our photometric sample using the equivalent 
width of the NB921 fiux; (total) ~276 h^^^ Mq yr-^ Given the survey volume, 217200 

Hq^ Mpc^, we nominally obtain a star-formation rate density (SFRD) of psFR.upper — 1-3 x 10"^ 
/io.7 Mq yr~^ Mpc"'^. However, this value should be regarded being as an upper bound obtained 
from our study. Since we estimated the fraction of LAEs in our photometric sample, /(LAE) 
~ 0.45, it seems reasonable to assume that 45% of the LAE candidates are at z ^ 6.6. We then 
obtained the corrected total SFR, (total) ^ 276 x 0.45 ~ 124 /i(^? Mq yr~\ This 

gives psFR — 5.7 x 10""^ /iq.t Mq yr^^ Mpc~^, being consistent with our previous study (Kodaira 
et al. 2003). Recently, Kurk et al. (2004) made a grism survey for LAEs at 2 ^ 6.4 - 6.6 and 
then identified a LAE at z = 6.518. Their new type of survey also gave a similar value of SFRD: 
psPR ^ 5.0 X 10-4 ho.7 Mq yr-i Mpc'^. 

We were also able to obtain a lower bound of the total SFR using the observed Lya fiuxes 
obtained with our spectroscopy. Using the data given in table 5, we obtained 5'Fi?*^P®'^ (total) 
~ 49.4 Hq^ Mq yr"^ This gives a lower bound of SFRD, psFR,iowcr ^ 2.3 x 10"^ Hq.y Mq yr"^ 
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Mpc"^. 

It is interesting to compare this value with previous estimates between z = and 2; ~ 6. In 
figure 16, we compare this star-formation rate density with those of previous studies compiled 
by Ajiki et al (2003) under the same cosmological parameters as those adopted in this paper. 
In addition, we also show the results obtained from LAE surveys by Cowie and Hu (1998), 
Kudritzki et al. (2000), Ouchi et al. (2003), and Ajiki et al. (2003) and from LBG surveys by 
Giavalisco et al. (2004b) and Dickinson et al. (2004). 

Our SFRD together with those based on the LAE surveys tend to give smaller values 
with respect to those obtained from LBG ones at 2; ~ 3 - 6. It should be reminded here 
again that we applied neither any reddening correction nor integration by assuming a certain 
luminosity function for the LAEs. Previous studies have shown that the SFR based on the Lya 
luminosity is smaller by a factor of two or more than that based on the UV continuum (Hu 
et al. 2002; Kodaira et al. 2003). Kodaira et al. (2003) also showed that the SFR based on 
the rest-frame optical continuum gives a value several times as high as that based on the Lya 
luminosity. More recently, Reddy and Steidel (2004) have shown from an analysis of GOODS 
multiwavelength data (see for GOODS, Giavalisco et al. 2004a) that the SFR based on X-ray 
data is higher by a factor of 5 than that based on the UV luminosity. In addition, integration 
using a luminosity function of LAEs also increases the SFR by a factor of two or more (see 
Ajiki et al. 2003). Although we do not still have a well-defined Lya luminosity function based 
on a large sample of LAEs (e.g., Santos et al. 2004; Hu et al. 2004), the SFR derived here for 
our sample of LAEs at z ^ 6.6 may be smaller by a factor of several than the real value; i.e., 
PsFR is on the order of 10~^ Hqj Mq yr~^ Mpc~^, or higher. 

Finally, we comment on SFRDs based on so-called i'-dropout samples. In addition to the 
follow-up optical spectroscopy of LAE candidates at z ^ 6.6 (this work), we also conducted op- 
tical spectroscopy of a sample of z'-dropouts found in the SDF (details will be given elsewhere). 
During this follow-up spectroscopy, we found a LAE at z = 6.33 (Nagao et al. 2004). Its z'-band 
flux is dominated by strong Lya emission, rather than the stellar UV continuum. If such strong 
LAEs could be involved in a photometric sample of i'-dropouts, we would overestimate SFRD. 
This issue shall be addressed in our forthcoming papers. 
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Table 2. Photometric properties of the LAEs . 
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(27.14) 


25.77 


l^ZD.U { ) 


1 9n 
izu 


Q7 


QT^T^ T1Q9/I9^ /l-i-979/(in 
oUr J ioz4z0.4+z / z4iU 


\9Q R9 


(27.58) 


25.57 


•^Ofi AO 
>"Z0.4Z 


\i 00 
> iyy 


00 


QT^T? Tl '?9/l9^ Q-l-97/l'?9/1 
our J ioZ4Z0.yn"Z / 4oZ4 


^98 R9 


>27.81 


25.56 


^Of\ AO 
P'ZD.4Z 


^971 
^z < 1 


Oi) 


our J loZ4Z0.U-t-Z / oDUO 


('98 

t^Zo. OD j 


(27.43) 


25.76 


^ZD.4Z 


^909 

^zyz 


AO 


STIT^ Tl '?94'?4 '?-l-974n^fi 

our J XOZ404.0^Z i 4UOU 


^98 fi9 


>27.81 


25.71 


^9(^ 49 


^0x0 


A 1 


'nT^T? Tl '?9/I "l^ n-U97'?Q^7 
our J loZ4oO.U-!-Z / oyO ( 


^98 R9 


>27.81 


25.45 


^OP AO 
P*ZD.4Z 


^ooy 


AO 
^Z 


CinTT Tl '?9/l "IR ^-^-97999*? 
our J ioZ4oD.o-rZ / zzzo 


^98 R9 


>27.81 


25.90 


^ZD.4Z 


^990 

^zzy 


/I Q 


oUr J ioz44U. z+z / zooo 


>Z(5.0z 


>27.81 


25.83 


>z0.4z 


>Z(54 


44 


SDF J132443.4+272633 


>28.62 


26.88 


25.37 


(25.89) 


163 


45 


SDF J132444.4+273942 


>28.62 


(27.29) 


25.41 


>26.42 


>299 


46 


SDF J132445.6+273033 


(28.20) 


26.97 


25.20 


(26.16) 


276 


47 


SDF J132447.7+271106 


>28.62 


26.92 


25.37 


(25.95) 


174 


48 


SDF J132449.5+274237 


>28.62 


(27.46) 


25.88 


>26.42 


>171 


49 


SDF J132450. 7+272160 


>28.62 


>27.81 


25.60 


>26.42 


>289 
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Table 2. (Continued.) 



■^0 




^9S fi9 


(97 7f\^ 


9^ P,7 




^9^9 


Ol 




^98 fi9 


(97 A^^^ 

t .^u } 


9^ 7I^ 


^9fi 49 




oz 


oJJr J ioZ4oo.O+z /oyio 


>ZO.DZ 


(Z ( .DU J 


ZO.OD 


>ZD.4Z 


>Ot)U 


53 


SDF J132458.0+272349 


>28.62 


(27.08) 


24.71 


>26.42 


>273 


54 


SDF J132503.4+273838 


(28.23) 


>27.81 


25.94 


>26.42 


>207 


55 


SDF J132506.4+274047 


>28.62 


>27.81 


25.95 


>26.42 


>193 


56 


SDF J132516. 7+272236 


>28.62 


(27.54) 


25.36 


>26.42 


>341 


57 


SDF J132528.0+271328 


>28.62 


(27.76) 


25.71 


>26.42 


>242 


58 


SDF J132533.4+271420 


>28.62 


>27.81 


25.65 


>26.42 


>336 



The LAEs are not detected in B, V, and R; see, for 3ct upper limits in these bands, table 1. 
'''The sky position, q;(J2000) and (5(J2000), is given in the name. 

■tAB magnitude in a 2" diameter. The magnitudes between the 1 a and 2 cr detection levels are put in parentheses. 
§LAE identified in Kodaira et al. (2003). 



Table 3. Journal of optical spectroscopy. 



MOS ID 


Setting* 


Exposure Time 


seeing 


Object IDt 






(s) 


(") 




Mask 1 


300R+O58 


19800 


0.5-0.8 


5, 7 


Mask 2 


300R+O58 


10800 


0.7-0.8 


4, 6 


Mask 4 


Echelle+z' 


14400 


0.6-0.7 


1,2,3 


Mask 5 


300R+O58 


12600 


0.6-1.0 


8 


Mask 6 


Echelle+z' 


13500 


0.5-0.7 


9 



Grating + order-cut filter. 

^Object ID numbers that correspond to the ones given in tables 2, 4, 5, 6 and 7. 
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Table 4. Spectral identification of emission-line objects. 



ID 


/red //blue 


log /red/ /blue 


1 


3.60±2.08 


0.556±0.251 


2 


1.22±0.43 


0.086±0.153 


3 


1.48±0.32 


0.170±0.094 


4 


1.33±0.41 


0.124±0.134 


5 


1.29±0.40 


0.111±0.135 


6 


1.33±0.18 


0.124±0.059 


7 


1.41±0.46 


0.149±0.142 


8 


1.15±0.54 


0.061±0.204 


9 


1 78±n 43 


25n±0 1 05 


12 


0.63±0.47 


-0.2G1±0.324 


16 


0.47±0.62 


-0.328±0.573 


17 


0.86±0.44 


-0.151±0.222 


18 


0.86±0.78 


-0.151±0.394 


19 


0.54±0.57 


-0.268±0.459 



Table 5. Star-formation properties of the nine LAEs. 



No. 




Name 


rP-(Lya) 






FWHM 








(lO^i® erg cm^^) 


(lO'*^ /ig 2 erg s-i) 


{h-^ Mq yr-i) 


(A) 


(km s 


1 


SDF 


J132352. 74-271622 


7.31 


3.56 


3.24 


14.0 


458 


2 


SDF 


J132353.1-I-271631 


8.22 


4.00 


3.64 


10.6 


347 


3 


SDF 


J132408.3+271543 


16.60 


8.13 


7.39 


7.5 


245 


4 


SDF 


J132415. 7-^273058 


19.60 


9.55 


8.69 


10.2 


334 


5 


SDF 


J132418. 3-^271455 


10.95 


5.41 


4.92 


< 9.0 


< 293 


6 


SDF 


J132418.4-I-273345 


19.43 


9.35 


8.51 


7.8 


256 


7 


SDF 


J132432. 5+271647 


6.17 


3.05 


2.77 


9.7 


316 


8 


SDF 


J132518. 8+273043 


7.29 


3.60 


3.28 


6.5 


212 


9 


SDF 


J132522. 3+273520 


15.48 


7.69 


7.00 


5.5 


179 
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Table 6. Full width at half maximum size of the nine LAEs. 



No. 




Name 


FWHM 


FWHAIcor* 


FWHMcor 








(") 


(") 


(K^ kpc) 


1 


SDF 


J132352. 7+271622 


0.81+0.04 


0.39+0.04 


2.1+0.2 


2 


SDF 


J132353. 1+271631 


0.95+0.03 


0.63+0.03 


3.5+0.1 


3 


SDF 


J132408.3+271543 


0.92+ 0.02 


0.59+0.02 


3.2+0.1 


4 


SDF 


J132415. 7+273058 


1.02+0.02 


0.73+0.02 


4.0+0.1 


5 


SDF 


J132418.3+271455 


0.84+0.02 


0.45+0.02 


2.5+0.1 


6 


SDF 


J132418.4+273345 


0.84+0.04 


0.45+0.04 


2.5+0.2 


7 


SDF 


J132432.5+271647 


0.93+0.04 


0.60+0.04 


3.3+0.2 


8 


SDF 


J132518.8+273043 


1.02+0.03 


0.73+0.03 


4.0+0.2 


9 


SDF 


J132522.3+273520 


0.94+0.03 


0.62+0.03 


3.4+0.2 



*Deconvolved with a PSF size of 0."71. 
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Table 7. Star-formation properties of our LAE candidates. 



No. 




Name 










S'Fi?(UV) 


1 


SUl' 


J132352.7H 


^271622 


4.4 


<6.4 


4.0 


<9.U 


2 


SDF J132353.H 


-271631 


4.6 


13.9 


4.2 


19.4 


3 


SDF J132408.3+271543 


9.2 


25.0 


8.4 


35.0 


4 


SDF J132415. 7+273058 


13.7 


30.1 


12.4 


42.2 


5 


SDF J132418.3+271455 


9.4 


<6.4 


8.6 


<9.0 


6 


SDF J132418.4+273345 


6.6 


<6.4 


6.0 


<9.0 


7 


SDF J132432.5H 


L 971 


4.8 


<6.4 


4.4 


<9.0 


8 


SDF J132518.8H 


:~z / oU^o 


4.4 


11.2 


4.0 


15.6 


9 


SDF J132522.3H 




7.1 


13.9 


6.5 


19.5 


10 


SDl' 


J 132406. 5H 


-271634 


3.4 


13.9 


3.1 


19.4 


11 


SDF 


J132410.8H 


-271928 


4.3 


21.0 


3.9 


29.4 


12 


SDF 


J132417.9H 


r271746 


5.4 


<6.4 


4.9 


<9.0 


13 


SDF 


J132428.7H 


^273049 


4.6 


<6.4 


4.2 


<9.0 


14 


SDF 


J132500.9H 


-272030 


4.1 


<6.4 


3.7 


<9.0 


15 


SDF 


J132515.5H 


-273714 


5.5 


<6.4 


5.0 


<9.0 


16 


SDF 


J132518.4H 


-272122 


15.4 


8.4 


14.0 


11.7 


17 


SDF 


J132519.9H 


-273704 


5.1 


<6.4 


4.7 


<9.0 


18 


SDF 


J132520.4H 


h273459 


4.0 


10.5 


3.7 


14.6 


19 


SDF 


J132521.H 


-272712 


5.1 


12.9 


4.6 


18.1 


20 


SDF 


J132525.3H 


-271932 


2.4 


14.6 


2.2 


20.4 


21 


SDl' 


J132338.4H 


-274652 


4.5 


<6.4 


4.1 


<9.0 


22 


SDF 


J132338.6H 


h272940 


4.1 


<6.4 


3.8 


<9.0 


23 


SDF 


J132342.2H 


h272644 


5.4 


<6.4 


4.9 


<9.0 


24 


SDF 


J132343.2H 


r272452 


10.9 


7.5 


9.9 


10.4 


25 


SDF 


J132347.7+272360 


4.5 


<6.4 


4.1 


<9.0 


26 


SDF 


J132348.9H 


h271530 


5.0 


<6.4 


4.5 


<9.0 


27 


SDF 


J132349.2+273211 


3.6 


<6.4 


3.3 


<9.0 


28 


SDF 


J132349.2+274546 


6.9 


<6.4 


6.3 


<9.0 


29 


SDF 


J132353.4+272602 


6.7 


11.3 


6.1 


15.9 


30 


SDF 


J132357.H 


h272448 


3.3 


8.0 


3.0 


11.2 


31 


SDF 


J132401.5H 


-273837 


3.6 


<6.4 


3.3 


<9.0 


32 


SDF 


J132402.6H 


-274653 


3.2 


13.9 


2.9 


19.5 


33 


SDF 


J132410.5H 


-272811 


4.0 


<6.4 


3.7 


<9.0 


34 


SDF 


J132419.3H 


-274125 


5.8 


<6.4 


5.3 


<9.0 


35 


SDF 


J132422.6H 


r274459 


5.2 


9.2 


4.7 


12.9 


36 


SDF 


J132424.2H 


h272649 


2.6 


8.8 


2.4 


12.3 


37 


SDF 


J132425.4+272410 


3.2 


<6.4 


2.9 


<9.0 


38 


SDF 


J132425.9H 


h274324 


4.4 


<6.4 


4.0 


<9.0 


39 


SDF 


J132425.0H 


h273606 


4.7 


<6.4 


4.3 


<9.0 


40 




J132434.3H 


h274056 


5.1 


<6.4 


4.6 


<9.0 


41 


SDF 


J132435.0H 


h273957 


5.4 


<6.4 


5.0 


<9.0 


42 


SDF 


J132436.6H 


h272223 


3.7 


<6.4 


3.3 


<9.0 


43 


SDF 


J132440.2H 


h272553 


4.6 


<6.4 


4.1 


<9.0 


44 


SDF 


J132443.4H 


h272633 


4.2 


10.4 


3.9 


14.5 


45 


SDF 


J132444.4H 


r273942 


4.8 


<6.4 


4.4 


<9.0 


46 


SDF 


J132445.6H 


h273033 


5.6 


8.1 


5.1 


11.3 


47 


SDF 


J132447.7H 


h271106 


4.3 


9.9 


3.9 


13.8 


48 


SDF 


J132449.5H 


h274237 


2.7 


<6.4 


2.5 


<9.0 


49 


SDF 


J132450.7H 


h272160 


4.6 


<6.4 


4.2 


<9.0 


50 


SDF 


J132455.4H 


^271314 


4.0 


<6.4 


3.7 


<9.0 


51 


SDF 


J132455.8H 


-274015 


3.2 


<6.4 


2.9 


<9.0 


52 


SDF 


J132458.5H 


h273913 


10.4 


<6.4 


9.5 


<9.0 


53 


SDF 


J132458.0H 


r272349 


4.4 


<6.4 


4.0 


<9.0 


54 


SDF 


J132503.4H 


h273838 


3.3 


<6.4 


3.0 


<9.0 


55 


SDF 


J132506.4H 


r274047 


3.1 


<6.4 


2.8 


<9.0 


56 


SDF 


J132516. 7+272236 


5.5 


<6.4 


5.0 


<9.0 


57 


SDF 


J132528.0H 


-271328 


3.9 


<6.4 


3.5 


<9.0 


58 


SDF 


J132533.4H 


-271420 


5.4 


<6.4 


4.9 


<9.0 



'The UV luminosity at A = 1260 AT 
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Fig. 1. The z' — NB921 color for a galaxy without a Lya emission line shown as a function of the redshift. 
The inset shows a close up for a redshift range between z = 6.0 and z =- 7.0. 
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Fig. 2. Color-magnitude diagram between z' — NB921 and A^i?921. All objects detected down to the 
apparent magnitude of NB921 = 26.5 in the -/Vi3921-selected catalog are shown. The horizontal red 
solid line corresponds to the color of z' — NB921 = 1.0 and the vertical red one corresponds to the So- 
limiting magnitude of A^i?921 = 26.0. The dashed curve shows the distribution of Str error in the color of 
z' ~NB921. The 58 LAE candidates in the photometric sample arc shown by filled circles, the twenty ones 
in the spectroscopic sample by filled circles with a large open circle, while the nine LAEs confirmed by our 
spectroscopy (see subsection 3.1) are shown with the same symbol but with a red large open circle. Data 
points without detection at z' are shown with a upper arrow; note that one a z' magnitude is adopted. 
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-1 t ^^^x^^^ u , ^ 

-10 12 3 

ijzz: 

Fig. 3. Diagram between z' — NB921 and i' — z' . The meanings of the symbols are the same as those in 
figure 2. See the text for the color evolution of a model galaxy (blue curve). 
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Fig. 4. Distributions of A^_B921 magnitudes for the photometric sample (58 objects; dotted histograms), 
the spectroscopic sample (20 objects; dashed histograms), and the spectroscopically-confirmed LAE sample 
(9 objects; solid histograms). 
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Fig. 5. Spatial distribution of the 58 objects given in our photometric sample (see table 2). The 20 LAE 
candidates in the spectroscopic sample are shown with a black open circle. The nine confirmed LAEs are 
shown with a red large open circle. 
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Fig. 6. Frequency distribution of log /rcd//biuo for the sample of 9 LAEs (solid-line) and 5 single-line 
objects (dotted- line) . 
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Fig. 7- 1. Observed optical spectra of nine LAE candidates (No. 1 SDF J132352. 7+271622). 
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Fig. 7- 2. Observed optical spectra of nine LAE candidates (No. 2 SDF J132353.1+271631). 
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Fig. 7- 3. Observed optical spectra of nine LAE candidates (No. 3 SDF J132408.3+271543). 
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Fig. 7- 4. Observed optical spectra of nine LAE candidates (No. 4 SDF J132415. 7+273058). This object 
was already reported in Kodaira et al. (2003). 
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Fig. 7- 5. Observed optical spectra of nine LAE candidates (No. 5 SDF J132418. 3+271455). This object 
was already reported in Kodaira et al. (2003). 



32 




Fig. 7- 6. Observed optical spectra of nine LAE candidates (No. 6 SDF J132418.4+273345). 
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Fig. 7- 7. Observed optical spectra of nine LAE candidates (No. 7 SDF J132432. 5+271647). 



34 




91 9150 9Z0a 



Wavelength (A) 



Fig. 7- 8. Observed optical spectra of nine LAE candidates (No. 8 SDF J132518.8+273043). 
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Fig. 7- 9. Observed optical spectra of nine LAE candidates (No. 9 SDF J132522.3+273520). 
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Fig. 8- 1. Observed optical spectra of five single-line objects (No. 12 SDF J132417.9+271746). 
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Fig. 8- 2. Observed optical spectra of five single-line objects (No. 16 SDF J132518.4+272122). 
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Fig. 8- 3. Observed optical spectra of five single-line objects (No. 17 SDF J132519. 9+273704). 
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Fig. 8- 4. Observed optical spectra of five single-fine objects (No. 18 SDF J132520.4+273459). 
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Fig. 8- 5. Observed optical spectra of five single-line objects (No. 19 SDF J132521. 1+272712). 
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Fig. 9. Thumbnail images of the four LAE candidates. North is up and east is left. The field of view is 
16" X 16", and the diameter of the circle is 8". 
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Fig. 10. High-resolution (PSF = 0."71) NB921 images of the nine LAEs. The number given in each 
panel corresponds to that given in tables 2, 4, 5, 6, and 7. 
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Fig. 11. Surface brightness radial distribution of tlie nine LAEs. The point-spread function is shown by 
a dashed blue curve and the best fit profile by black one in each panel. The number given in each panel 
corresponds to that given in tables 2, 4, 5, 6, and 7. 
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Fig. 12. Histograms of L™'^s''(Lya) for the photometric sample (58 LAE candidates) and for the nine 
LAEs. The latter is shown by the thick line. The bins correspond to Icr {— 1.4 x 10^^ /iq y erg s~^), 2(T, 
3a, and 4cr, of U'^'^^'^{Lya); see table 7 in more details. 
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Fig. 13. Comparisons between /^^'^'^{Lya) and /™'^8°(LyQ;) for the nine LAEs. 
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Fig. 14. Histograms of _L^(UV) at A = 1260 A for the pliotomctric sample (tlie 58 LAE candidates) and 
for tlie nine LAEs. The latter is shown by the thick line. The bins correspond to la (= 0.64 x 10^^ /iq 7 
erg s~^ A~^), 2(7, and 3a of Ly(UV); see table 7 in more details. 
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Fig. 15. Comparisons between S'Fi?™'^s<=(LyQ;) and SFR{\JV) for tiie 58 LAE candidates. The nine 
LAEs arc shown by filled circles. The inset shows a close up to see crowded data points for S F R™^'^^" (Lya) 
V.7 ^^O yr"^ and for 5Fi?(UV) = 5 15 /i^^ Mq yr'^. 
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Fig. 16. Star-formation rate density shown as a function of the redshift. Our data points are shown 
with big symbols; psFR = 5.5 x 10^"' ho,j Mq yr~^ Mpc~'^ is shown by the filled circle. Pgpp''^ is shown 
by an upper half-filled circle while Pg°^^' by an upper half-filled one. The other data sources are Gallego 
et al. (1995 - plus), Tresse and Maddox (1998 - star), Fujita et al. (2003b - open square), Treyer et al. 
(1998 - filled square), Lilly et al. (1996 - open circles), Connolly et al. (1997 - open triangles), Madau et 
al. (1996 - filled triangles), Steidel et al. (1999 - open inverse triangles), Iwata et al. (2003 - filled inverse 
triangle), Giavalisco et al. (2004b - open diamonds), Bouwens et al. (2004 - filled diamond), Stanway et 
al. (2004b - open hexagon), and Ouchi et al. (2004 - filled hexagons). For reference, we also show the 
results of previous Lya searches at z ~ 3 - 6 by Kudritzki et al. (2000 - double square), Cowie and Hu 
(1998 - double circle), Fujita et al. (2003a - double triangle), Ouchi et al. (2003 - double diamond), and 
Ajiki et al. (2003 - double hexagon). 
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